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Abstract 

Three applied-field MPD thruster geometries were tested with argon propellant to 
establish the influence of electrode geomety on thruster performance. The thrust increased 
approximately linearly with anode radius, while the discharge and electrode fall voltages 
increased quadratically with anode radius. All these parameters increased linearly with 
applied-field strength. Thrust efficiency, on the other hand, was not significantly 
influenced by changes in geometry over the operating range studied, though both thrust and 
thermal efficiencies increased monotonically with applied field strength. The best 
performance, 1820 sec I S]) at 20% efficiency, was obtained with the largest radius anode at 
the highest discharge current (1500 amps) and applied field strength (0.4 Tesla). 
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applied magnetic field strength, T 
electron diffusion coefficient 
electron charge, C 
electric field 

electron current density. A/cm 2 
discharge current, A 
Boltzmann's constant, J/K 

Coulomb logarithm 
mass flow rate, g/s 

electron mass, kg, and number density, cm -3 
anode and cathode radii, cm 
electron Larmor radius, cm 
thrust, N 

electron temperature, eV 
Alfven critical ionization velocity, m/s 
discharge and electrode fall voltages, V 
linear distance 

Debye length, cm 

thrust and thermal efficiencies 

electron number flux 

electron mobility and permeability of free space 
electron cyclotron frequency, s* 1 
electron - ion collision frequency, s* 1 
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Introduction 


Analyses of lunar and Mars missions utilizing nuclear electric propulsion systems 
have shown that the thruster specific impulse should be between 5000 and 8000 seconds 
with an efficiency over 50 percent 1 . In addition, the propulsion system must deliver a total 
impulse of approximately 1 x 10 9 N-s while operating at average power levels between 2 
and 10 MW*A While MPD thrusters may be able to fulfill these requirements, the large 
disparities between the mission requirements and the capabilities of current technology have 
forced a multi-faceted approach to address these issues, including direct measurements in 
quasi-steady thrusters at multi-megawatt power levels and scaled testing of steady-state 
thrusters at lower powers. This paper presents preliminary results of one aspect of the 
research directed at establishing the feasibility and practicality of MPD thruster technology. 

The highest efficiencies and specific impulses are reported for multi-megawatt 
quasi-steady MPD thrusters with applied magnetic fields using hydrogen propellant, with 
reported values reaching 70% at 7000 seconds I sp 3 . Quasi-steady self-field thruster 
efficiency is generally lower 3 ' 5 . Thrusters utilizing other propellants (Ar, He, NH 3 , N 2 , 
and Ne) have still lower performance, with peak values of approximately 30 - 40% and 
specific impulses up to 4000 seconds 2 . Results for steady-state thrusters operated at 
powers up to 600 kW cover a similarly wide range, with the highest performance reported 
with an applied field thruster operated on lithium at powers below 30 kW (7000 sec I sp , 
60% efficiency 6 ). Self-field thruster testing at powers up to 600 kW has been done with 
argon 7 , with a maximum performance of less than 25% efficiency at 1500 seconds I sp . 

Results to date have been obtained with a broad range of thruster geometries, with 
power levels ranging from tens to thousands of kilowatts, propellant flow rates from 0.005 
to 10 g/s, and a variety of propellants. The very range of the MPD thruster research makes 
it at present difficult to conclusively establish the performance trends for either self or 
applied-field thrusters. This problem is compounded by the lack of detailed theoretical 
models to guide the design of more efficient thrusters. Previous modelling efforts for both 
self-field 8 ' 9 and applied-field 16 ' 12 thrusters have shown substantial promise, but they have 
not been verified or implemented for a variety of thruster conditions to indicate directions 
for research. For example, the importance of the shape and strength of the magnetic field 
or the electrode size and geometry are not presently understood. The combined lack of an 
extensive database of directly measured performance and detailed models has precluded 
significant advances in MPD thruster performance. 

This paper presents preliminary results of an effort to systematically study the 
impact of geometry and operating condition on MPD thruster performance. Following a 
description of the test facility and the thruster geometries, a summary of the performance 
measurements will be given, including measurements of power lost to the electrodes. 
Simple models are then used to interpret the observed trends. 


Experimental Apparatus 


Test Facility 

The MPD thrusters were mounted on a thrust stand in a 3 m diameter, 3 m long test 
section of a 7.6 m diameter, 21 m long main vacuum chamber. The 3 m diameter gate 
valve between the test section and the main tank was opened after inserting the thrust stand 
into the test section. The test stand is shown in Fig. 1 . The tank pressure was maintained 
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below 0.07 Pa during testing at argon flow rates up to 0.14 g/s. Thruster power was 
provided by a series-parallel network of 66 kW welding supplies with a maximum 
capability of 3000 A at 88 V. The propellant for these tests was argon. Propellant flow 
rates and thrust levels were calibrated to within 5% and 2%, respectively. In-situ thrust 
calibrations were done using a motor driven weight suspension apparatus. Details of the 
thrust stand and its calibration, tares, etc., are given in Ref. 13. 

The thruster discharge was initiated using a combination of a 1 kV, 2 amp power 
supply connected across the thruster electrodes and a short duration gas pulse to 
temporarily raise the thruster chamber pressure. The main power supplies were protected 
from the high voltage pulse using blocking diodes in series with the supplies. Generally a 
1-2 second gas pulse followed by a 0.1 second voltage pulse was sufficient to ignite the 
discharge. A short current overshoot , lasting 1 - 3 seconds, preceded stabilization at the 
preset current level. 

Water cooling for the thruster and applied-field magnet was provided by two 
closed-loop heat exchangers. The water pumps supplied up to 0.8 1/sec at a pressure of 
IMPa. Water flow rates were measured using turbine flow meters calibrated to within 2%. 
Further descriptions of the facility are found in References 13-15 . 

MPD thrusters and Applied-Field Magnet 

The MPD thruster geometries tested in this work are shown in Fig. 2. These 
coaxial thrusters consisted of a central 1.27 cm diameter, 2 % thoriated tungsten cathode 
with an annular copper anode. The anodes were water cooled along their entire length, 
while the cathodes were water cooled behind the thruster backplate. Two cylindrical 
anodes and one flared anode were tested in this work. To permit comparisons between 
geometries the "effective anode radius" of the flared anode was defined as the average of 
the upstream and exit plane radii. While this appears arbitrary given our limited knowledge 
of the current distribution, reasons for the choice will become clear in the next sections. 
Thruster dimensions are given in Table 1. Geometries A and B were modified versions of 
commercially available plasma sprayers. These designs used O-ring seals, which failed at 
high input powers. These rubber seals were eliminated in geometry C, which extended the 
thruster power range. 

The applied-fieid magnet consisted of a solenoid made from 46 turns of water 
cooled copper tubing (1 .27 cm O.D, 0.89 mm wall ) in a 14 cm ID, 15 cm long spool. The 
magnet had four layers of ten coils each, and a final layer with six coils. Thrusters were 
placed in the magnet such that the thruster exit plane was aligned with the magnet face to 
reduce the heat transfer to the solenoid from the thruster plume. The magnetic field at the 
magnet center reached a maximum of 0.4 Tesla with 1500 amps through the solenoid. A 
typical plot of axial applied-field strength vs. axial position along the centerline is given in 
Fig. 3. The axial field strength at the thruster exit plane is plotted vs. radius in Fig. 4. In 
this paper, the axial field strength at the centerline of the thruster exit plane is used as the 
reference value. 


Experimental Results 

The discharge voltage and thrust for the three thruster geometries, operated at a 
discharge current of 1000 amps and an argon flow rate of 0.14 g/s, are plotted versus the 
applied magnetic field strength (at the centerline of the thruster exit plane) in Figures 5 and 
6, respectively. The discharge voltage increased approximately linearly with applied-field 
strength for all geometries. Linear least-squares fits to these data showed that for the two 
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cylindrical geometries the slope of the discharge voltage versus applied field strength curve 
increased quadratically with the anode radius, while the flared anode has a lower slope than 
that implied by this relationship. The curve-fit results are summarized in Table 2. For 
example, the slope for geometry A, with a radius of 1.3 cm, was 37.2 V/Tesla, while for 
geometry C, with a radius twice that of geometry A, this slope was 149.8 V/Tesla, or four 
times the slope observed with geometry A. This scaling relationship will be discussed 
further in the next section. No clear correlations were found between the discharge voltage 
slope and the anode or cathode length, or with the discharge current. 

As shown in Fig. 6, the thrust level also increased approximately linearly with 
applied-field strength for all geometries, and the slope is again strongly dependent on anode 
radius. As seen from the values in Table 3, however, the quantitative relationship between 
the slope and the anode radius is not as apparent. To further study this geometric effect, 
two tests were conducted with the cylindrical geometries (A and C) in which both the 
discharge current and the applied-field strength were varied. As shown in Fig. 7, the thrust 
increased linearly with the product of the discharge current and the applied-field strength 
(JdxB) for both geometries. In addition, the data in Figure 7 show that the slope of the 
thrust- JdxB curve is directly proportional to the anode radius. 

The thrust measurements resulted in the specific impulse values, defined by 


shown in Fig. 8. The maximum value of 1820 seconds was reached with the largest 
diameter anode at the highest discharge current and magnetic field strength. 

The thrust efficiencies, defined by 


T 2 

2mV d J d 


( 2 ) 


are shown in Fig. 9 for a discharge current of 1000 A and argon flow rate of 0.14 g/s. 
The efficiencies for all three geometries were comparable, though at low applied-field 
strengths the smallest anode had a significantly higher efficiency. The efficiencies 
approach 20 % for all geometries at high field strengths. It is clear that the efficiency of 
the larger anode (geometry C) depended more strongly on field strength than the smaller 
anodes (geometries A and B). In addition, the leveling off in the efficiency vs. applied 
field curve observed for the smallest anode (geometry A) was not seen with the larger 
anodes. 


The principal cause of the low thrust efficiencies was electrode power deposition. 
The fraction of the total power deposited into the electrodes was established using 
calorimetric measurements of the water cooling the electrodes. These measurements were 
used to establish the thermal efficiency, defined by 


^ = 1 


Electrode Power 

vS 


(3) 


which is the fraction of the input power going into the plasma. The water flow rate to the 
electrode cooling passages was adjusted to maintain a minimum of six degrees Celsius 
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temperature change between inlet and outlet. This ensured electrode power measurements 
to better than 10% with a thermocouple accuracy of 0.56 degrees Celsius. 

Results of these measurements are shown in Fig. 10, where the thermal efficiency 
is plotted vs. the applied field strength. The thermal efficiencies for the two smaller radius 
thrusters are only weakly dependent on the magnetic field strength, and there is substantial 
scatter in the data for geometry A. Geometry C, with the largest anode radius, showed the 
strongest dependence on applied field strength, with values increasing from 25% at 0.1 T 
applied-field to 40% at 0.4 T. It is clear that for all geometries over 50% of the input 
power was deposited into the electrodes. 

The influence ol the magnetic field on electrode power deposition is more clearly 
shown using an "effective electrode fall voltage", defined as the total power deposited into 
the electrodes divided by the discharge current. Results for the three geometries, shown in 
Figure 11, indicate a linear rise in the effective electrode fall with magnetic field. The 
smallest geometry has both the lowest magnitude electrode fall and the slowest rise with 
applied-field strength. The highest value, 44 volts, is reached with the largest thruster 
(geometry C) at the highest field strength. Linear least-squares fits for the two cylindrical 
geometries, summarized in Table 4, indicate that the slopes increased quadratically with the 
anode radius. Once again, using the effective radius for the flared anode thruster yields a 
slope too high for this relationship. 


Discussion 


The above data have several implications. First, since thrust is approximately 
proportional to the anode radius, and the discharge voltage is proportional to radius 
squared, the thrust efficiency should be approximately independent of anode radius, as 
observed. However, the thermal efficiencies for the three geometries were significantly 
different for low applied-field strengths, with the larger anode thruster losing more power 
to the electrodes. This indicates that the larger anodes must recover more of the energy 
deposited into the plasma. Thus, while the overall power balance does not appear to 
change, the relative magnitudes of the power sinks do. 

The measured values of specific impulse, while low when compared to mission 
requirements, are much higher than expected for simple self-field acceleration. Using the 
standard self-field thrust equation with an effective anode to cathode radius ratio ^ 


T = 


4n 


In 


V 

V r c J 


eff 


(4) 


for the geometries tested yields specific impulses of 52, 75, and 102 seconds, respectively, 
for the test conditions. These values are orders-of-magnitude below those measured. 
While at low applied-field strengths thermal effects may be important with the small anode 
geometries, the high specific impulses observed for high applied-fields cannot be explained 
by thermal acceleration of argon. It is clear that the dominant plasma acceleration was due 
to the applied-field. 

The theoretical model that most clearly lends itself to interpreting the above 
measurements is that of Fradkin'H He assumed that the applied-field acceleration was 
dominant, that rotational motion was converted to axial motion with constant efficiency, 
and that the plasma followed a magnetic field flux tube, and solved the equations of 
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continuity and conservation of angular momentum and energy to establish the parametric 
behavior of applied-field MPD thrusters. Fradkin experimented with only one thruster 
geometry, and so was not able to verify the model predictions with regards to anode shape. 
While highly simplified, some of the general trends predicted by the model agree with 
those observed in this work. Specifically, Fradkin predicted that the electromagnetic thrust 
and discharge voltage scale as 


and 


T =W BV - 


1 — 


U) 

k t ») 


v d =v c+ ^i 

d c 2m 


(5) 

( 6 ) 


where V c is the voltage due to ohmic heating and the electrode falls. The thrust scaling 
arises from the dependence of the applied-field torque on the product of the field strength 
and the discharge current. Note that for the geometries studied, the second term in the 
brackets is small. The linear increase of thrust with anode radius is due to the linear 
decrease in current density with anode radius and the quadratic increase in cross-sectional 
area. The voltage scaling arises from the back-emf (UxB) term in Ohm’s Law: the 
velocity is proportional to BJ(jr a (from thrust), so the product UxB scales as B 2 Jdr a . The 
second power of the radius arises from the linear increase in path length with anode radius. 

The linear increase in thrust with BJdr a matches the results shown in Figure 7. 
However, while the discharge voltage plotted in Figure 5 does increase quadratically with 
anode radius, it does not appear to increase quadratically with the magnetic field strength 
and testing revealed that it is not sensitive to the discharge current. This same discrepancy 
was observed by Fradkin, who attempted to explain it by invoking the Alfven critical 
velocity hypothesis. For a velocity limited case, the voltage should scale as 


V d = V c +U c Br, 


r 2 ) 

r J 


(7) 


where U c is the Alfven critical velocity. This has the same behavior as our measurements 
except for the radius dependence. Thus, while the model correctly predicts thrust scaling 
behavior, there is a problem with the voltage prediction which must be resolved by future 
models. 


The discrepancy between the geometry scaling for the cylindrical anodes and that 
for the flared anode is most probably due to a non-uniform current distribution in the latter 
geometry. To substantiate this conclusion, the "true effective anode radius" for the flared 
anode (geometry B) was calculated assuming the scaling laws of V ~ r a 2 , T ~ r a and V e ~ 

r a 2 , and using the measured slopes to calculate r a . The results were uniformly smaller than 
the geometric average radius by between two and three millimeters, indicating that the 
current attached preferentially on the upstream side of the anode in the flared anode 
thruster. 


The linear increase in the effective electrode fall voltage with applied magnetic field 
indicates that the sheath is magnetized. This occurs whenever the electron gyro radius is 
comparable to or smaller than the Debye length, and the electron cyclotron frequency is 
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much larger than the electron- ion collision frequency. These ratios are calculated from 17 


and 



7.5xl0~ n 

eB 


V m e n e 


co M eBT c 15 

v e _j 2xlO _6 m e n e ln(A) 


( 8 ) 

(9) 


Note that the first relation is independent of temperature. Values for argon at several 
magnetic fields and electron densities are given in Table 5. The Coulomb logarithm and the 
ion charge were given constant values of ten and one (singly ionized), respectively, for 
these estimates, so the results should be within a factor of two of actual values. The 
electron densities bracket measurements made using electrostatic probes close to the 
thruster exhaust plane 18 . It is clear that for most of the operating conditions studied here 
the electron Larmor radius is comparable to the Debye length and that electrons cycle about 
magnetic field lines many times before undergoing a collision. These estimates substantiate 
the conclusion that the magnetic field impedes electron motion in the sheath, and they have 
important implications for methods used to model electrode processes in applied-field 
thrusters. 

The linear rise in effective fall voltage with applied-field strength can be explained 
in part by considering the diffusion of electrons into the anode. Previous studies of steady- 
state thrusters have concluded that the anode fall is the dominant source of electrode power 
deposition, so that the data presented above are most likely representative of processes 
occuring near the anode 19-21 . The electron flux into the anode can be written as 22 

r e = -D e ^-^ e n c E (10) 

dy 


Using the Bohm diffusion coefficient 17 for D e and rearranging using Einstein's relation 
between the diffusion coefficient and the electron mobility, we obtain 


E = - 


1 

en e 


16j e B + kT e 


dy J 


(ID 


which shows that, to a f irst order, the electric field required to maintain a constant electron 
current increases linearly with magnetic field strength. Note that classical diffusion would 
predict a quadratic increase in electric field with magnetic field strength, a dependence 
which is clearly not observed in these measurements. Limitations of this argument include 
its neglect of ion current and our lack of knowledge of how the other terms vary with 
magnetic field strength. 


Conclusions 


Prelimary studies of the influence of applied-field thruster geometry on performance 
showed that the anode radius had a large effect on thrust, discharge voltage and electrode 
fall voltage for a given operating point. Specifically, thrust appeared to increased linearly 
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with anode radius, and the discharge and electrode fall voltages increased quadradically 
with anode radius. These quantities increased linearly with applied field strength over the 
operating range covered in this study, and thrust was proportional to the product of the 
applied field strength and the discharge current. These effects are such that the thrust 
efficiency is not significantly influenced by thruster geometry, though both thermal and 
thrust efficiency increased monotonically with magnetic field strength. The results show 
that while thruster geometry did not dramatically affect the overall thruster power balance, 
the magnitudes of the dominant power sinks changed with geometry. It is important to 
note that the operating range studied here did not include regions where self-field effects 
contribute significantly to plasma acceleration. 
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Tahle 5: Ratio of Larmor radius to Debve length and gyro frequency 
to electron-ion c ollision frequency 
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Figure 2: MPD Thruster Geometries 
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AXIAL DISTANCE FROM MAGNET CENTER, CM 

Figure 3: Axial applied-field strength vs. axial position along 

centerline. 



RADIUS, CM 


Figure 4: Axial applied-field strength at thruster exit plane vs. 

radius. 
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MAGNETIC FIELD AT EXIT PLANE, T 

Figure 5: Discharge voltage vs. magnetic field strength at 

thruster exit plane centerline (Jd=1000A, m=0.14g/s). 



MAGNETIC FIELD AT EXIT PLANE, T 
Figure 6: Thrust vs. applied magnetic field strength at exit 

plane centerline (Jd=1000A, m=0.14g/s). 
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SPECIFIC IMPULSE, SEC 



DISCHARGE CURRENT x MAGNETIC FIELD, A-T 

Figure 7: Thrust vs. the product of the dischage current 
and the magnetic field strength (Jd=1000A, m=0.14g/s). 



DISCHARGE CURRENT x MAGNETIC FIELD, A-T 
Figure 8: Specific Impulse vs. the product of the discharge 
current and the magnetic field strength (Jd=1000A, m=0.14g/s). 
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THERMAL EFFICIENCY 
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MAGNETIC FIELD AT EXIT PLANE, T 

Figure 9: Thrust efficiency vs. magnetic field strength at the 

thruster exit plane centerline (Jd = 1000 A, m = 0.14g/s). 
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MAGNETIC FIELD AT EXIT PLANE, T 

Figure 10: Thermal efficiency vs. magnetic field strength at the 

thruster exit plane centerline (Jd=1000 A, m = 0.14 g/s). 
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MAGNETIC FIELD AT EXIT PLANE, T 

Figure 11: Effective electrode fall vs. magnetic field strength 

at the thruster exit plane centerline (Jd=1000 A, m = 0.14 g/s). 



18 




IVI/NSA Report Doc 

Nalional Aeronautics and I IC^UI L UUk 

Space Administration 

1. Report No. NASA TM-103259 2 Governmen 

AIAA-90-2669 
4. Title and Subtitle 

Geometric Effects in Applied-Field MPD Thrusters 


Report Documentation Page 


2. Government Accession No. 


3. Recipient’s Catalog No 


5. Report Date 


6. Performing Organization Code 


7. Author(s) 

R.M. Myers, M. Mantenieks, and J. Sovey 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


8. Performing Organization Report No. 

E-5692 

10. Work Unit No. 

506-42-31 

1 1 . Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Prepared for 21st International Electric Propulsion Conference sponsored by AIAA/DGLR/JSASS, Orlando, 
Florida, July 18-20, 1990. R.M. Myers, Sverdrup Technology, Inc., Lewis Research Center Group, 2001 
Aerospace Parkway, Brook Park, Ohio 4-4142; M. Mantenieks and J. Sovey, Lewis Research Center. 

16. Abstract 

Three applied-field MPD thruster geometries were tested with argon propellant to establish the influence of electrode 
geometry on thruster performance. The thrust increased approximately linearly with anode radius, while the 
discharge and electrode fall voltages increased quadratically with anode radius. All these parameters increased 
linearly with applied-field strength. Thrust efficiency, on the other hand, was not significantly influenced by 
changes in geometry over the operating range studied, though both thrust and thermal efficiencies increased 
monotonically with applied field strength. The best performance, 1820 sec I sp at 20% efficiency, was obtained 
with the largest radius anode at the highest discharge current (1500 amps) and applied field strength (0.4 Tesla). 


17. Key Words (Suggested by Author(s)) 

Electric propulsion 
Magnetoplasmadynamic thrusters 
Plasma accelerators 


18. Distribution Statement 

Unclassified — Unlimited 
Subject Category 20 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21 . No. of pages 
20 


NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 




